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Scaling, Finite Size Effects, and Crossovers of the Resistivity and Current- Voltage 
Characteristics in Two-Dimensional Superconductors 
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(Dated: March 1, 2013) 

We revisit the scaling properties of the resistivity and the current- voltage characteristics at and 
below the Berezinskii-Kosterlitz-Thouless transition, both in zero and nonzero magnetic field. The 
scaling properties are derived by integrating the renormalization group flow equations up to a scale 
where they can be reliably matched to simple analytic expressions. The vortex fugacity turns out to 
be dangerously irrelevant for these quantities below Tc, thereby altering the scaling behavior. We 
derive the possible crossover effects as the current, magnetic field or system size is varied, and find 
a strong multiplicative logarithmic correction near Tc, all which is necessary to account for when 
interpreting experiments and simulation data. Our analysis clarifies a longstanding discrepancy 
between the finite size dependence found in many simulations and the current-voltage characteristics 
of experiments. We further show that the logarithmic correction can be avoided by approaching the 
transition in a magnetic field, thereby simplifying the scaling analysis. We confirm our results by 
large scale numerical simulations, and calculate the dynamic critical exponents z, for relaxational 
Langevin dynamics and for resistively and capacitively shunted Josephson junction dynamics. 

PACS numbers: 74.40.-n,74.78.-w,64.60.Ht 



Fluctuation effects can be very strong in low- 
dimensional systems and may radically alter the mean 
field picture of phase transitions. A well known ex- 
ample is that of two-dimensional (2D) superfluids or 
superconductors, where phase fluctuations of the com- 
plex order parameter ip = ipge^^ destroy long range or- 
der at all nonzero temperatures. Despite this, a super- 
fluid/superconducting phase with algebraic order, finite 
superfluid stiffness, and zero resistivity, still exists at low 
temperature. This is separated from the high tempera- 
ture disordered phase by a transition - the Berezinskii- 
Kosterlitz-Thouless (BKT) transition - caused by the 
thermal unbinding of vortex-antivortex pairs [lH3| . The 
properties of 2D superconductors have been studied in- 
tensely in recent years [J, |5| and continue to receive much 
interest due to the relevance for cuprate superconduc- 
tors with their layered structure. Furthermore, advances 
in fabrication enables studies of single or few atomic 
layer thick superconductors, which offer great potential 
for precise tests against theories and simulations [6|. In 
this paper we explore the possible scaling behaviors and 
crossover effects that may occur as function of current, 
magnetic field, and system size. These results are con- 
firmed by numerical simulations and used for an accurate 
determination of the dynamic critical exponent for two 
different equations of motion. 

Transport measurements are perhaps the best way to 
experimentally study the properties of 2D superconduc- 
tors. One of the hallmarks of the BKT transition is the 
nonlinear current-voltage (IV) characteristics E ~ J"*^^^ 
at and below Tc, with a temperature dependent expo- 
nent a{T) [71, y. The exact form of the temperature 
dependence of the exponent a{T) has been subject to 
some debate [9[, but there is a quite general agreement 
that it should take on the value a(TBKT) = 3 right at 



the transition. According to the conventionally accepted 
theory developed by Ambegaokar, Halperin, Nelson and 
Siggia (AHNS), a{T) = cahns = 1 + 27rJ^(r)/2T, where 
J{T) = h^ps{T)/2m is the superfluid stiffness and PsiT) 
the (fully renormalized) superfluid areal density [7|, |8|. 
This result has been contested by Minnhagen et al. 
(MWJO) [3 who arrived at the alternative expression 
a(r) = aMwjo = 2TrJ'{T)/T—l using scaling arguments. 
Both yield a = 3 at Tbkt — '^J{T)/2. Alternatively one 
may try to describe the data using a Fisher-Fisher-Huse 
(FFH) scaling formula [l^ 



E = J^-^-^-^f (J^'^-i/T), 



(1) 



where £{■) is a scaling function and ^ the correlation 
length. This leads also to a power-law, but leaves a = 
z -|- 1 as a free fitting parameter related to the dynamic 
critical exponent z (d = 2 is the dimension). In 2D, 
however, fits of experimental data to Eq. ([Ij easily give 
surprisingly large values a > 6 ll[, although more rea- 
sonable values ~ 3 have also been obtained 15|. This, 
however, highlights the difficulty in using Eq. ([T]) with- 
out additional assumptions. In any case it remains chal- 
lenging to decide which of the scenarios described above 
is correct based only on experiments. One may instead 
resort to computer simulations to try to settle the con- 
troverse, but this too is no easy task. Usually, simulation 
data is analyzed using finite size scaling formulas based 
on Eq. (HI), with the diverging correlation length ^ cut 
off by the system size L, yieldin g E ~ JL^~°' for small 
J. Most [12-16J (but not all [H Si) simulation studies 
appear to favor the value aMWJO- Interestingly, Refs. Il8l, 
obtain agreement with both the AHNS and MWJO ex- 
pressions in different regimes and for different boundary 
conditions. At the same time, the validity of the FFH 
scaling formula Eq. ^ is still an open question, as is the 



scaling behavior in the presence of an apphed magnetic 
field. 

The main contribution to the scaling behavior of the 
resistivity and IV characteristics conies from the free vor- 
tex density np of unbound vortex pairs. These can be 
either thermally excited or induced by an applied mag- 
netic field or a current. The BKT transition occurs when 
thermally excited vortex-antivortex pairs unbind. Since 
only the motion of free vortices dissipate energy, the re- 
sistivity should be proportional to the free vortex density 



where $o is the flux quantum and fi^ 
Bardeen-Stephen vortex mobility. 

Conventionally, the free vortex density np = Tif 
is calculated from a rate equation [7[ 



(2) 
27rf^p„/$g is the 



drip 
dt 



= r - An 



F"'F^ 



(3) 



where F = Az^e"'^"'*/"^ is the pair generation rate and A 
the recombination rate. Here z = e~^=/^ is the vortex fu- 
gacity, and Ef. ^^ J the vortex core energy. The potential 
barrier to overcome in order to create a pair of free vor- 
tices has two terms, one which depends logarithmically 
on their separation r, and one with a linear dependence 
due to the applied current C/cff(r) w 27rJ'ln(r/ao) — J$o^, 
where oq w ^o is a short distance cutoff of the order of 
the Ginzburg-Landau coherence length. (From now on 
we set flo — 1.) Optimizing gives r* sa 2ttJ /^qJ and 
U^s = t7eff(r*) w -2ttJ [ln( J$oao/27rj7) + 1]. The sta- 
tionary solution to Eq. ([3]) gives 



np ~ 2ze 



-U^iillT 



2zJ 



2-nJ 12T 



(4) 



and, with E = p{J)J, the result a = cahns- 

There are several ways in which the above picture may 
need to be modifled. First, interactions between vor- 
tices except those constituting the pair are completely ne- 
glected. Screening of the vortex interaction from bound 
vortex-antivortex pairs can be taken into account by us- 
ing the fully renormalized value of the stiffness J7(T) in 
place of the bare one. In a finite system the vortices may 
enter and exit the system at the boundaries and Eq. ^ 
will acquire more terms describing these processes. Ac- 
counting for a realistic geometry and nonuniform cur- 
rent distribution can lead to a rather complicated behav- 
ior [19[. In simulations one usually avoids surface effects 
by using periodic boundary conditions (PEG). Finite size 
effects, however, become visible when r* = 2ttJ'/^oJ > 
L, leading to a crossover to ohmic behavior at low cur- 
rents, with a characteristic size dependent resistivity. 
Another issue is that the rate equation ([3]) presumes that 
density fluctuations are small, which is true for large 
systems, but not for small enough systems with area 



L^ < l/np- In the latter regime the constraint of vortex- 
antivortex neutrality (enforced when using PEG [20]) in- 
stead leads to F/A = {n^rip) « ('^_f) +L^'^ {''^f)^ which 
is dominated by the second term, i.e.. 



The same expression follows from a low fugacity expan- 
sion of the neutral Goulomb gas, which only involves even 
powers of z. Also note that an applied perpendicular 
magnetic field B will lead to a net density of free vor- 
tices An = nj — rip = B/^q, such that 



An + AripUp 



An' 



- Az'e- 



-U^tt/T 



(6) 



A more systematic approach to take into account in- 
teraction effects, is to first integrate the renormalization 
group (RG) flow up to the scale where one of the coupling 
constants becomes large of 0(1) and only then match the 
theory to simple approximate expressions similar to the 
ones discussed above. The RG flow equations are most 
easily expressed in the Goulomb gas language using the 
rescaled temperature and fugacity variables, a: = 1 — ^, 
y = 211 z. To lowest order in x and y they read [3|, [21 1 



dx „ 9 

-n = 2y', 



dy „ 



(7) 



where ^ = ln6 is the logarithm of the scale factor h. The 
resulting RG flow obeys x' — y' = C' , where 



\C\ = J\xl-yl\^c^W^^T\ 



(8) 



is a constant determined by the initial conditions. Ee- 
low Tc we have C^ > and the RG flow ends up on a 
critical line x = — C < 0, y = as £ — > cxi. Above Tc, 
C' < and the flow will eventually diverge to -t-oo. The 
BKT transition occurs at T — Tc, where the flow follows 
the separatrix x = —y. In order to describe the vari- 
ous crossovers we need the explicit solutions [2l|, y{£) = 
C/sinh(2C(£-£o)) for T < T^ y{i) = l/{2i-2io) for 
T = Tc, and j/(£) = -\C\/sin{2\C\{i ~ to)) for T > T^. 
In terms of 6 = e we have 



yib) = 
y(h) 



-2C 



2Cib/bo) 
1 



21n(V&o)' 



{T < Tc), (9) 

{T = Tc), (10) 



where 6o — e^" is fixed by the initial conditions. Near 
Tc, where \C\ < j/q, we have to a good approximation 
c2 « Ayo/nJ, T, « ^^/2(1 + yo), Iq ~ -l/2?/o- Further 
below Tc, where C ^ yo, we have instead C ~ xq, so that 

yib) ^ yob-^"" , (r«T,). (11) 

Note also that C = -x{b ^ oo) = ttJr{T)/2T - 1 is 
directly related to the fully renormalized superfluid stiff- 
ness Jr{T). 



The free vortex density, being the vortex density which 
remains after the eUmination of all bound pairs, is only 
rescaled by the RG transformation and therefore has scal- 
ing dimension 2, i.e, np ^ b~^. As a function of system 
size L, magnetic flux density B, current J, x, y, and 
possibly other perturbations it therefore transforms as 



nF{xo,yo,L,B,J,...) = 

b-^npixib), y{b),Lb-\Bb^, Jb, . 



(12) 



under the RG. A similar equation holds for the resistivity 
Eq. ^. Most theories assume that the vortices undergo 
ordinary diffusion However, we are not aware of any ar- 
gument which prevents the renormalization of the vortex 
mobility fiy in Eq. ^. Hence, we allow for an anomalous 
dimension /j,^ ~ &^~^, with a dynamic critical exponent z 
not necessarily fixed to 2, such that the resistivity trans- 
forms as 

p{xo,yo,L,B,J,...) = 

b-'p{x{b),y{b), Lb-\Bb^, Jb, . . .). (13) 

An FFH scaling formula follows from Eq. P^ if p flows 
smoothly to a nonzero constant as & — >■ oo. This is the 
case above Tc, where the flow must be stopped at a scale 
when X ~ y '^ 0(1), yielding the Debye-Hiickel expres- 
sion np ~ l/27r^^, where ^+ ~ exp(7r/2cV^ — ^c) is 
the correlation length above Tc. This is, however, not 
the case in zero magnetic field at and below Tc, where 
y — 2ttz — >■ 0, because np vanishes in this limit. In other 
words, the fugacity is dangerously irrelevant for np and p 
in this case. Instead the right hand side of Eqs. (|12p -(fT3 |) 
must be matched to one of Eqs. (H])-®. At the matching 
scale b the barrier Uoff in Eq. (jl]) or ((5)) has reduced to 
zero, and we are left with three different possibilities: In 
zero magnetic field np{b) ^ y{b) or y^{b) depending on 
boundary conditions and system size, while for nonzero 
field np{b) « y/B^b*/^l + y^ib)/?:^. This wih turn out 
to have profound consequences for the scaling of many 
quantities. 

We first discuss the finite size scaling of the linear re- 
sistivity in zero magnetic field. The RG flow must then 
be stopped a.t b = L. Under the RG all length scales, in- 
cluding the system size, shrink by a factor b so that the 
effective system size becomes L' ^ L/b = 1. The system 
must therefore be matched to Eq. ([5]) when using periodic 
boundary conditions, or to ([4]) when using open bound- 
ary conditions. For PBC we thus get p{L) ^ L~^y'^{L), 
and by using Eqs. (P)-([TT|). the hmiting cases 

( ]^-z+4-27,Jn/T (L> f ) 

\L Vln {L/bQ), (£ < ^-), 



where ^_ sa exp(l/2C) « exp{l/2c^/T~^) is the cor- 
relation length below Tc, defined as the scale on which 
x{b) has approximately reached its asymptotic value — C. 
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FIG. 1. (Color online) Langevin dynamics, (a) Resistivity p 
vs system size L at different temperatures in a magnetic field 
B — $o/L^ ■ The dotted lines are guides for the eyes, and the 
full green curve at T = Tc is a x^'ht (using L = 16-120) to 
the power-law p '^ L~^, giving z — 2.22. (b) As in (a), but 
for zero magnetic field. The full green curve at T = Tc is a 
X^-fit (using L = 16-80) to p ~ L''^-^^ /(InL - eof with fixed 
z = 2.22, giving io ~ —2.71. Insets: The effective exponent 
Zcff VS temperature, obtained from power-law fits. Note how 
ZcH is almost constant below Tc in (a). 



The power-law appearing in this expression agrees with 
the finite size scaling of MWJO [1^ if one assumes 
z = 2. On the other hand, for open boundary condi- 
tions p{L) ^ L^^y{L), or 



P{L) 



j^-z+2-ttJr/T 

L-y\niL/bo), 



(i<e-), 



(15) 



which, for z = 2, would be consistent with the AHNS 
scaling. The finite size scaling at Tc, where ^_ = cxd, has 
in both cases, strong multiplicative logarithmic correc- 
tions. 

The situation in a nonzero magnetic field is different. 
The magnetic field is a relevant perturbation, which de- 
stroys superconductivity by introducing a finite density 
of free vortices even at low temperature. We can, how- 
ever, still approach the transition by scaling down the 
magnetic field with the system size, holding BL^ = N^q, 
the net number of flux quanta, fixed. (This is easy in a 
simulation, but more difficult in an experiment.) Con- 
sider, e.g., the case A^ = 1. Stopping the RG fiow at 
b - L 
L- 



y^^o/B and matching to Eq. ([6]) then gives 
The leading scaling behavior 



p^ L Vl + 2/'WA'- 
thus remains a temperature independent power-law with 
exponent z in contrast to the zero field case (with a weak 
additive correction decreasing with system size). 

The finite size scaling formulas derived above are well- 
suited for the analysis of numerical simulations. We have 
performed simulations of the 2D XY model, defined by 
the Hamiltonian H = —JJ2(ij) cos{9i — 6j — Aij), using 
two types of dynamics, relaxational Langevin dynamics 
and resistively and capacitively shunted Josephson junc- 
tion (RCSJ) dynamics (in the overdamped limit) [22l |. 
The resistivity was calculated from the equilibrium volt- 
age fluctuations using a Kubo formula, with a sampling 



time of 10^-10^ time units per datapoint. For an accu- 
rate determination of z we apply a weak magnetic field 
B = ^q/L^ so that the system contains exactly one vor- 
tex irrespective of system size. This minimizes the in- 
fluence of the logarithmic correction near Tc, allowing 
us to fit the data for T < Tc to the simple scaling law 
p{L) ~ L~^. We plot, in Fig. HJa), p vs L calculated 
using Langevin dynamics on a log-log scale for a range 
of temperatures including T^ {T^ « 0.892J (H). The 
data at and below T^ do indeed follow a power-law with 
a temperature independent exponent z w 2.22 ± 0.05. 
In contrast, the zero field data shown in Fig. [Ijb) fol- 
low different power-laws at different temperatures. Right 
at Tc the data is very well fitted by Eq. ()14|) with z 
fixed to 2.22. The value of 4 = In^o ~ -2.7 ob- 
tained by the fit compares well with the theoretical es- 
timate io « —1/22/0 ~ —2 obtained using the XY value 
2/0 = 27re^^<=/^, with Ec ~ Tr^J/2. Without knowing 
about the logarithmic correction one would fit the data 
at Tc to a pure power-law and draw the wrong conclu- 
sion. For our data this would give an effective exponent 
ZeS ~ 2.54, appreciably different from the true z. 

Figure [D shows similar plots for RCSJ dynamics. The 
resistivity for a system with exactly one vortex again fol- 
lows a power-law, but this time with z = 1.77 ± 0.05 at 
Tc- In zero field the data is well fitted to (fT4|) using the 
same z, with £o ~ —1.33 again in rough agreement with 
expectations, whereas a pure power-law fit would give a 
too large exponent z^s ~ 2.2. 

The values z « 2.22 and z « 1.77 for Langevin and 
RCSJ dynamics, respectively, are close to, but signifi- 
cantly different from the conventional value 2, and cor- 
respond either to subdiffusivc (z > 2) or superdiffusive 
{z < 2) vortex motion. 

The scaling behavior below Tc differs considerably in 
zero and nonzero magnetic field. As seen in the insets 
of Figs. [T] and [5] the resistivity with B = ^q/L"^ follows 
a power-law with practically temperature independent 
exponents in stark contrast to the zero field case. Previ- 
ous finite size scaling studies of p{L) (or E{J,L) in the 
ohmic regime) in zero field have obtained a temperature 
dependent power-law exponent below Tc in good agree- 
ment with the MWJO prediction 13l4l6l Il8 |. which is 
not surprising given (fT4)) and the smallness of z — 2. 

In a large or infinite system at zero magnetic field the 
RG flow must be stopped at a scale dictated by the ap- 
plied current, i.e., when J5 ~ Jo = 2ttJ'/^o. At this 
scale the matching condition is np ^ y and the nonlinear 
resistivity p{J) ^ J^y{b ~ Jo/ J) obtains from Eqs. ^- 
pT|) . We have the limiting cases 



.M^f 



jz+7tJr(T)/T-2^ 

JVln(Jo/J6o), 



Jo/ J > e- 
Jo/J < e- 



(16) 
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FIG. 2. (Color online) Overdamped RCSJ dynamics (Stewart- 
McCumber /3c = 2-n-IcR^C/^o = 0.25). (a) Resistivity p vs 
system size L at different temperatures in a magnetic field 
B = $o/i^. The full green curve at T = Tc is a x^-fit (using 
L — 16-120) to the power-law p ~ L~^ , giving z — 1.77. (b) 
As in (a), but for zero magnetic field. The full green curve at 
T^Tcisax^-fit (using L = 16-80) to p r^ L'^'^'' /{lnL-£o)' 
with fixed z — 1.77, giving £o — —1.33. Insets: As in Fig. [T] 



Close to Tc we find a strong multiplicative logarithmic 
correction. The crossover to the finite size induced ohmic 
behavior in Eq. p^ or ([15]) happens when JL < Jo. In 
addition one expects a high-current crossover to an ohmic 
regime when J ^ Jo- 
in the PBC case it is also possible to have an interme- 
diate regime where the matching is still done at a scale 
b « Jo/ J, but the effective system size is small enough 
that nF{h){L/hY ^ 1, so that np ~ y^. This would give 



piJ,L) ^ L^J 



2 jz-2+2ttJr/T 



J 



<L< 



ttJr/2T 

• (17) 



The power-law behavior at low currents below Tc is in 
agreement with the AHNS value if one assumes z — 2. 



Such an intermediate scaling regime was previously pro- 
posed in Ref. [181, using an entirely different approach. 

To summarize, we have obtained a coherent picture of 
the scaling behavior and crossover effects of the (non- 
linear) resistivity near and below the BKT transition, 
Eqs. p^ - pT|) . The finite size results depend sensitively 
on the boundary conditions and on whether a magnetic 
field is present or not. In the limit of large systems the 
IV exponent agrees with the AHNS result, with the mod- 
ification that we allow for the possibility that z ^2. For 
PBC, on the other hand, the finite size scaling agrees 
with MWJO. Our simulations suggest that z differs from 
2 and moreover that Langevin and RCSJ dynamics be- 
long to different dynamic universality classes 2J]. From 
a practical point we found it important to take into ac- 
count the logarithmic correction near Tc when analyzing 
finite size data. The same should hold true for experi- 
mental finite current data. Note, however, that to make 
quantitative comparisons with experiments it may be im- 
portant to consider effects of inhomogeneity and pinning, 
and to make realistic estimates of the temperature depen- 
dence of the bare parameters J7, y, e.g., using Ginzburg- 
Landau theory [2511. Finally, it should be noted that the 



only assumptions needed in our analysis is the low fugac- 

ity behavior of the zero magnetic field resistivity p ^ y 

or y^. It is highly likely that other quantities may be "\ 

affected in similar ways. 

We thank M. Wallin for comments on the manuscript. 
This work was supported by the Swedish Research rg- 

Council (VR) through grant no. 621-2007-5138 and the 
Swedish National Infrastructure for Computing (SNIC [10 
001-10-155) via PDC. 

[11 



* Ianan02@kth.sel 
^ jlidmar@kth.se 

[1] V. S. Berezinskii, Sov. Phys. JETP 32, 493 (1971)34, 610 
(1972) 

[2] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 
(1973) 

[3] P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987) 

[4] J. M. Repaci, C. Kwon, Q. Li, X. Jiang, T. Venkates- 
san, R. E. Glover, C. J. Lobb, and R. S. Newrock, 
,Phys. Rev. B 54, R9674 (1996)S. T. Herbert, Y. Jun, 
R. S. Newrock, C. J. Lobb, K. Ravindran, H.-K. Shin, 
D. B. Mast, and S. Elhamri, 57, 1154 (Jan 1998)L. Miu, 
D. Miu, G. Jakob, and H. Adrian, 73, 224526 (Jun 
2006)M. Mondal, S. Kumar, M. Chand, A. Kamlapure, 
G. Saraswat, G. Seibold, L. Benfatto, and P. Raychaud- 
huri, |Phys. Rev. Lett.| 107, 217003 (Nov 2011)T. I. Ba- 
turina, S. V. Postolova, A. Y. Mironov, A. Glatz, M. R. 
Baklanov, and V. M. Vinokur, EPL (Europhysics Let- 
ters) 97, 17012 (2012) 

[5] D. R. Strachan, C. J. Lobb, and R. S. Newrock, 
|Phys. Rev. B 67, 174517 (May 2003) 

[6] N. Reyren, S. Thiel, A. D. Caviglia, L. F. Kourk- 
outis, G. Hammerl, C. Richter, C. W. Schneider, 
T. Kopp, A.-S. Riietschi, D. Jaccard, M. Gabay, D. A. 
MuUer, J.-M. Triscone, and J. Mannhart, Science 317, 
1196 (2007)G. Logvenov, A. Gozar, and I. Bozovic, 
[Science (Ne w York, N.Y. ) 326, 699 (Oct. 2009), ISSN 
1095-9203J. T. Ye, S. Inoue, K. Kobayashi, Y. Kasahara, 



[12: 

[13; 

[14 
[15 
[16 
[17 

[is; 

[19 
[20 
[21 
[22; 



[23' 
[24; 

[25] 



H. T. Yuan, H. Shimotani, and Y. Iwasa, Nature mate- 
rials 9, 125 (Feb. 2010), ISSN 1476-1122 
V. Ambegaokar, B. I. Halperin, D. R. Nelson, and E. D. 
Siggia, Phys. Rev. Let t. 40, 783 (1978 )Phys. Rev. B| 21, 

1806 (1980) 

B. I. Halperin and D. R. Nelson, ^JLTP" 36, 599 (1979) 

P. Minnhagen, O. Westman, A. Jonsson, and P. Olsson, 

|Phys. Rev. LetE] 74, 3672 (1995) 

D. S. Fisher , M. P. A. Fisher, and D. A. Huse, 

|Phys. Rev. BJ 43, 130 (1991) 

S. W. Pierson, M. Friesen, S. M. Ammirata, J. C. Hun- 

nicutt, and L. A. Gorham, Phys. Rev. B 60, 1309 (Jul 

1999) 

K. Medvedyeva, B. J. Kim, and P. Minnhagen, 

|Phys. Rev. Bl 62, 14531 (2000) 

B. J. Kim, P. Minnhagen, and P. Olsson, |Phys. RevT^ 

59, 11506 (1999) 

L. M. Jensen, B. J. Kim, and P. Minnhagen, 

Phys. Rev. B 61, 15412 (2000) 

P. Minnhagen, B. J. Kim, and A. Gronlund, 

Phys. R ev. B, 69, 064515 (2004) 

H. Weber, M. Wallin, and H. J. Jensen, [Piys. Rev. B| 53, 

8566 (Apr. 1996) 

M. V. Simkin and J. M. Kosterlitz, |Phys. Rev.^ 55, 

11646 (1997) 

Q.-H. Chen, L.-H. Tang, and P. Tong, |Phys. Rev. Lett.| 

87, 067001 (2001)L.-H. Tang and Q.-H. Chen, 

Phys. Rev. B 67, 024508 (2003) 

A. Gurevich and V. M. Vinokur, |Phys. Rev. Lett.| 100, 
227007 (2008) 

By PBC we mean, here and in the following, any bound- 
ary condition which enforces vortex neutrality. 
J. M. Kosterlitz, Journal of Physics C: Solid State Physics 
7, 1046 (1974) 

Details of the simulation methods can be found in A. An- 
dersson and J. Lidmar, Phys. Rev. B 83, 174502 (May 

2011) 

P. Olsson, IPhys. Rev. B| 52, 4526 ( 1995) 

P. C. Hohenberg and B. I. Halperin, [Rev. Mod. Phys.| 49, 
435 (Jul 1977) 

L. Benfatto, C. Castellani, and T. Giamarchi, 
Ph ys. Rev. B| 80, 214506 (2009) 



